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ABSTRACT 

A comprehensive analysis of the available literature on 
the atmosphere of Jupiter is made, and two models are selected 
as being most likely. It is suggested that a simple probe to 
Jupiter might resolve the question of which atmosphere actually 
exists . 

The general ballistic entry equations are displayed, and 
specific solutions for direct and displaced orbital entries 
are made in two dimensions for Jupiter. It is found that a 
large dense vehicle, requiring the Saturn booster and utilizing 
a graphite ablative heat shield, can make a direct parabolic 
entry at angles up to 30° • It is doubtful that mid-course and 
terminal guidance exists which can give the' precise injection 
angles required by this mission. 

The feasibility of atmospheric determination by measure- 
ments of drag load factor, altitude, and heating rate is 
discussed. The general conclusion is that if a successful 
entry can be made, the atmosphere can be determined. The 
success of the entry depends on control of the entry angle. 
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OBJECT 



The object of this thesis is to examine the vehicle design 
parameters for an entry into the atmosphere of Jupiter, and to 
investigate the feasibility of atmospheric determination by a 
simple entry vehicle. 



CHAPTER 1 



THE ATMOSPHERE OF JUPITER 



1.1 Introduction 

The purpose of thrs thesis is to examine the design parameters 
for a ballistic entry vehicle into the atmosphere of Jupiter. The 
first step is to determine the most probable model of the Jovian 
atmosphere . 

The Jovian atmospheres are quite different from those of the 
Terrestrial planets. Jupiter is the closest and easiest to reach 
of these major planets, and therefore will provide Man's first 
experiment with such an atmosphere. 

The most generally accepted theory of planetary formation is 
Kuiper's protoplanet theory (ref l) . The Solar nebula was broken 
up by gravitational instability and turbulence into the several 
protoplanets and the Sun. Condensation was already present, as 
the temperature of the gas in this nebula was much lower than that 
of the present planets. The heavy solid particles ranging from 
icr-5 to io5 cm in diameter were drawn together to form the nucleus 
of the protoplanets. This only represented about one percent of 
the mass material. The heavier gases, notably oxygen, nitrogen, 
the noble gases, and the oxides of carbon, nitrogen and sulphur, 
were drawn toward the center and trapped. Many of the minerals 
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were carbonated, sulpha ted, and nitrated. 

The remaining 99 percent of the mass in the gaseous envelope 
proceeded to evaporate due to the temperatures generated by contrac- 
tion, and due to solar radiation on the nearer protoplanets. This 
envelope of Hydrogen, Helium, ammonia and methane completely evaporated 
from the Terrestrial planets. Some masses broke off due to the 
centrifugal force of high rotation, to later form the satellites of 
Earth and Mars. A few percent of the original envelope remained 
with the Jovian planets, however, and this accounts for their re- 
ducing hydrogen atmospheres. Jupiter retained approximately 10 per- 
cent of its protoplanet mass, and thus is the largest planet. 

Because of their higher velocities around the Sun, the outer planets 
lost more of their original envelope, thus they are smaller. It is 
believed that Pluto lost all of its gaseous envelope and therefore 
closely resembles Earth, except for its temperature. 

Further contraction of the Terrestrial planets forced percola- 
tion and exhalation of gases from within the hot core, and this 
accounts for the Earth’s oxidizing atmosphere. Such atmospheres as 
remain on the Terrestrial planets, then, are these secondary emissions 
of the heavy gases from within. Mercury lost most of its secondary 
atmosphere because of its size and proximity to the Sun. Venus 
apparently has lost most of its water due to solar proximity, and Mars 
has lost most of its atmosphere because of its small size. 

The evaporation from the major planets continued for a shorter 
time than it did from the Terrestrial planets. Kuiper estimates 10 
years for Earth, and 10'’ years for Jupiter. This is due to the lower 
temperature and greater distance from the Sun, and also to the larger 
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mass of Jupiter. The rate of evaooration depends on the log of 
the mass to radius ratio, and the temperature. Kuiper has shown 
that Jupiter lost considerable Hydrogen, and some Helium during 
the period of 10^ years when it could lose it. At the present time, 
Jupiter is losing a negligible amount of gas. We should note that 
the loss in H^ proceeds as the square of the loss of He, due to 
the masses of the two molecules. Therefore, contrary to the 
opinion of Urey (ref 2), considerably more Hy drogen could be 
lost than Helium. 

* 

We shall consider the atmosphere in more detail in the 
following sections. This should be our first experience with a 
reducing, light gas, "fuel" atmosphere. 

1.2 The Ability to Hold an Atmosphere 

According to the kinetic theory of gases (e.g., Jeans 
(ref 3) )» an atmosphere in approximate isothermal equilibrium 
has the following density profile: 

_ mgr / H \ 

P ' fe e* RT r + H (1-1) 



where r = radius of planet, H = height above datum, and 
= datum density. 

Jeans also states that for stability over astronomical 

Q 

periods (10 7 years), the RMS molecular velocity of the gas at 
the exosphere should be less than 20 percent of the escape 
velocity, V esc . 



V 

esc 



T 

Rj + H 



Ft. 

Sec. 



(1=2) 



For Jupiter, this works out to 64.6 km/sec. or 196,800 ft. /sec l 
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At the present time, with temperatures of the order of 86 °K, very few 
molecules at the escape level attain this velocity. However, during 
the protoplanet period, with temperatures of the order of 1000°K in 
the exosphere, and very low M/R ratios, efficient escape of would 
occur. According to Kuiper, the value of log M/R for Jupiter at the 
beginning of the evaporative period was -.39. The upper limit of 
log M/R permitting efficient escape is; +.18 for m = 2 (hydrogen), 
-.12 for m = 4 (helium), and -.77 for m = 18 (methane, ammonia). 

Thus we see that plenty of hydrogen, some helium, but no methane 
could have evaporated. 

The other Jovian planets have less ability to hold their 
envelopes, and it is estimated that Saturn holds 2.5 percent and 
Uranus and Neptune about 1.25 percent. On this basis alone, Uranus 
and Neptune would be expected to have little or no hydrogen, which 
is not the case. However, it is believed that the extreme cold 
temperatures on these planets may have slowed diffusion of H^ and 
He to the exosphere during the evaporation period. Their atmospheres 
also show a lack of turbulence which is present on Jupiter and 
Saturn. Their relative abundances cannot therefore be accurately 
predicted, but must be empirically determined. The usual assumption 
is that fractionation factors are some smooth function of m, the 
molecular weights of the constituents. 

1.3 Probable Constituents 

The composition of the protoplanets at the time of their 
formation is estimated to be in the solar or cosmic proportions. 
Accurate spectrographic analyses of the Sun and stars give the 
relative amounts of gaseous elements shown in Table I. 
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TABLE I 

Initial Composition of Gaseous Envelope 


Gas 


Amount 


Molecular Weight 


Weight 


i° Weight 


H 2 


6289 


2.016 


12,679 


63.5 


He 


1740 


4.003 


6,965 


34.9 


Ne 


5,9 


20.18 


119 


0.6 


H 2° 


3.75 


18.02 


68 


0.34 


nh 3 


3.0 


17.03 


51 


0.26 


ch 4 


1.4 


16.04 


22 


0.11 


A 


0.8 


36.3 


29 


0.15 


Totals 


8044 




19.933 


100 i 


Average molecular weight, m = 2.48 


r 


= 1.46 



Urey and others believe that the present remaining envelope 
of Jupiter contains hydrogen, helium and methane in these same 
proportions above the tropopause . To Kuiper it seems unlikely 
that the gases would remain so homogeneous, and would escape 
uniformly, leaving solar proportions intact, Kuiper holds that the 
gases would be fractioned in proportion to m^ , the molecular 
weight of each gas. These separate opinions lead to two different 
models of the atmosphere. To this writer, it seems more likely 
that Kuiper is correct, and this theory is supported by experi- 
mental evidence to be discussed later. However, we will carry 
both atmospheres through this thesis to provide a basis for 



6 



comparison of design parameters. 

Urey's atmosphere, model "a", is shown below in Table II. 



TABLE II 

Solar Proportion Atmosphere, model a 


Gas 

H 

2 

He 

CH> 

4 


Amount 
6 289 
1740 
1.4 


Molecular Weight 
2.016 
4.003 
16.04 


Weight 

12,679 

6,965 

22 


$ Weight 

64.4 

35.4 
0.2 


Totals 8030 




19,666 


100 i 


Average molecular weight, m = 2.47 


= 1.47 



Kuiper's atmosphere^ model b, is characterized by a reduction 
of Urey's model "a" by a factor such that the ratios of Hydrogen to 
Helium to Methane are: 

6289n 2 : 1740n : 1.4 

H 2 : He : CH 4 (l-3) 

The factor n here is the reduction due to evaporation, and as 
pointed out before, H 2 evaporates as the square of He. To make the 
present mass come out right, and with the assumption of about 10 $ 
of the protoplanet envelope remaining, n = .177* Therefore, the 
model b atmosphere is reproduced in Table III. 

We notice that to change the molecular weight ty even this 
small amount requires a considerable change in the composition. 
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TABLE III 

Reduced Proportion Atmosphere 


, model b 




Gas 


Amount 


Molecular Weight 


Weight 


# Weight 


h 2 


197 


2.016 


396 


23.9 


He 


308 


4.003 


1234 


74.7 


ch 4 


1.4 


16.04 


22.5 


1.4 


Totals 


5 06.4 




1652.5 


100.0# 


Average molecular weight, m = 3.26 


0 


= 1.56 



The pressures in the two models will vary considerably, as shown 
in the section on structure. These models are only for the gas 
above the tropopause, as it is assumed the heavier elements will 
have very low pressures at this height. Notice also that model b 
has He in excess of H 2 by 3 to 1, whereas in the solar proportions 
is over He ty 2 to 1 ty mass. Urey holds that such a great 
reduction in H 2 would be unlikely. Perhaps, then, the actual 
atmosphere is somewhere between these models, for surely some of 
the gas must have evaporated. We will carry both atmospheres 
through our calculations for comparison. 

1.4 Methods of Measurement 

The primary method used to determine atmospheric composition 
has been comparison of absorption spectrographs of the planet 
with emission spectrographs of the Sun. Absorption of a parti- 
cular frequency line in the Sun's spectrum indicates the presence 
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of a particular gas. The strength of the absorption indicates the 
amount. This method has several drawbacks, however. It is limited 
in scope because the planets are not primary emitters; only the light 
from the Sun can be re-radiated. Certain gases like nitrogen and 
helium have no absorbtion lines in the useful range of spectrography. 
The Earth's atmosphere shuts out most of the ultra violet wave- 
lengths. We also get interference from CO^ , H 2 0, and oxygen in 
the Earth's atmosphere. Most of these hindrances can be avoided 
with careful work, but the most difficult of all is to obtain the 
laboratory comparison spectra of the various gases for identification. 

Methane and ammonia bands were detected as early as 193? ty 
Wildt in the atmosphere of Jupiter. The ammonia band was much 
weaker than expected. It was then realized that the clouds on 
Jupiter were NH^ cirrus, and less than one cm NPT* of free NH^ 
existed above the cloudtops. According to Kuiper's protoplanet 
theory, the atmosphere should contain an abundance of H 2 and He, 
but no bands for these gases were identified. In 1938 Herzberg 
identified several lines around 7500 A° in the atmospheres of all 
the Jovian planets; and lines around 8270 A° in Uranus and Neptune. 

But it was not until 1952 that a laboratory experiment showed the 
8270 A° lines to be pressure-induced quadruple vibration-rotation 
modes of the H 2 molecule. To reproduce this on Earth required a 
path length of 80 meters, a temperature of 78°K gained by cooling 
with liquid nitrogen, and a pressure of 100 atmospheres. 

This illustrates the extremes in pressure, path length, and 
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cm NPT = height in cm of equivalent column of the gas on Earth 
at normal pressure and temperature. 



temperature that must be used to demonstrate on Earth the same 
lines found in nature's laboratory. It can be shown that the 
temperature on Jupiter at high pressures is too high to allow 
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this excitation. But can be confirmed by its presence on 
Uranus and Neptune. The lines at 750 0 A° are present on 
Jupiter, but have not to date been identified. 

A rare but extremely important method of constitutive 
determination was the occultation of a star by Jupiter in 
November 1952. The rate of diminution of light from the star as 
it passes behind the planet is a function of the apparent motion, 
temperature, surface gravity, and mean molecular weight of the 
gases in the atmosphere. Since the first three of these parameters 
are fairly well known, the mean molecular weight can be found by 
photometric measurement of the rate of light diminution. 

It was by this method that Baum and Code (ref 4) measured 
the m of Jupiter in 1952. The light curve best fit an atmosphere 
with m = 3*3 • This was important confirmation that the primary 
constituents were and He, and CH^ and NH^ as indicated spec- 
tographically. It further indicates an abundance by weight of 
He over ^ by 3^1 as predicted by Kuiperl This is why this 
author prefers Kuiper's model b atmosphere to Urey's model a. 

Other methods of atmospheric measurement include ultra- 
violet and infrared absorption by the atmosphere, polarization, 
and radio emission. Although radio emission in a band around 
22 me has been detected from Jupiter, no conclusive evidence of 
composition has been found. Astronomers believe that this 
radiation is "static", caused by vast turbulence below the 
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clouds on Jupiter. According to Peek (ref 5) this turbulence gains 
its energy from heat within the planet, and from the fast planetary 
rotation. The turbulence accounts for the meteorological formation 
of the bands on the planet similar to the prevailing winds on Earth. 
This will be discussed more fully under Structures. 

1.5 Structure of the Atmosphere 

The general structure of both models will conform to the 
following assumptions: 

1) A spherically symmetric atmosphere 

2) Isothermal equilibrium in the stratosphere 

3) Dry adiabatic lapse rate in troposphere 

4) Moist (NH^) adiabatic lapse rate below the cloud tops. 

From these assumptions, the thermal models can be constructed 

as shown in figure 1. The isothermal temperature given by Peek 
(ref 5) is 86°K. The equilibrium temperature for solid' NH^ from 
the phase diagram is l65°K to l68°K, depending on pressure. Zero 
conditions refer to the tops of the clouds. The measurements of 
Jupiter's radius have always been referred to the cloud tops, so 
we will use this as zero altitude. 

The dry adiabatic lapse rate is given by the expression: 

d-4) 

g ^ 

This yields for Jupiter: 

Ijp = 2.47)= 2.6 °K/km 

r (m = 3.26^= 4.0 °K/km 

With this information, the depth of the cloud tops below the 
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tropopause is easily calculated. On the way down a level is 
passed which is at the radiometric temperature of Jupiter (102°K) . 
This is the level which the thermocouple sees at 8-l4p. It is also 
the level where the pressure is approximately one atmosphere. The 
emissivity of the gases above this point is so low that we assume 
most of the energy radiated into space comes from below this 
level. Between here and the cloud tops there is less than one 
centimeter NPT of NH^, and radiation from below the clouds is 
effectively blocked. This explains the rather weak ammonia 
absorption spectrographs from Jupiter. 

The density ratio for the adiabatic portion of the atmos- 
phere follows; 

^ (1 - 5) 

' o O o 



The density, p^, is calculated from the perfect gas law, which 
should apply at these temperatures and pressures. 



/°o 




R T 
g o 



(1-6) 



The pressures at the cloud tops are fixed by the (H^ + He/CH^) 
abundance ratio used for model a and b„ 



P (m = 2.47) = 24 earth atmospheres 

P (m = 3.26) = 2 earth atmospheres 

Here we point out again that model b is more probable. 

The 8270 A° absorption due to quadruple excitation is absent 
from Jupiter. But if a pressure of 24 atm prevailed in the 
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clouds on Jupiter. According to Peek (ref 5) this turbulence gains 
its energy from heat within the planet, and from the fast planetary 
rotation. The turbulence accounts for the meteorological formation 
of the bands on the planet similar to the prevailing winds on Earth. 
This will be discussed more fully under Structures. 

1.5 Structure of the Atmosphere 

The general structure of both models will conform to the 
following assumptions: 

1) A spherically symmetric atmosphere 

2) Isothermal equilibrium in the stratosphere 

3) Dry adiabatic lapse rate in troposphere 

4) Moist (NH^) adiabatic lapse rate below the cloud tops. 

From these assumptions, the thermal models can be constructed 

as shown in figure 1. The isothermal temperature given by Peek 
(ref 5) is 86°K. The equilibrium tanperature for solid' NH^ frcm 
the phase diagram is l65°K to l68°K, depending on pressure. Zero 
conditions refer to the tops of the clouds. The measurements of 
Jupiter's radius have always been referred to the cloud tops, so 
we will use this as zero altitude. 

The dry adiabatic lapse rate is given by the expression: 

r=^ (*fr> (1-4) 

This yields for Jupiter: 

fa = 2A7j= 2.6 °K/km 

r^ m = 3.26/= 4.0 °K/km 

With this information, the depth of the cloud tops below the 
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tropopause is easily calculated. On the way down a level is 
passed which is at the radiometric temperature of Jupiter (102°K). 
This is the level which the thermocouple sees at 8-l4p, . It is also 
the level where the pressure is approximately one atmosphere. The 
emissivity of the gases above this point is so low that we assume 
most of the energy radiated into space comes from below this 
level. Between here and the cloud tops there is less than one 
centimeter NPT of NH^, and radiation from below the clouds is 
effectively blocked. This explains the rather weak ammonia 
absorption spectrographs from Jupiter. 

The density ratio for the adiabatic portion of the atmos- 
phere follows: 




The density, is calculated from the perfect gas law, which 
should apply at these temperatures and pressures. 



/o 




R T 
g o 



( 1 - 6 ) 



The pressures at the cloud tops are fixed ty the (H^ + He/CH^) 
abundance ratio used for model a and b. 



P (m = 2.47) = 24 earth atmospheres 

P (m = 3*26) = 2 earth atmospheres 

Here we point out again that model b is more probable. 

The 8270 A° absorption due to quadruple excitation is absent 
from Jupiter. But if a pressure of 24 atm prevailed in the 
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visible portion of the atmosphere, we probably could see this band. 

Above the tropopause isothermal conditions are assumed to 
prevail. Therefore the density follows the expression: 



Most references give the scale height of Jupiter as about 
60,000 feet, but such a figure is based on a mean temperature of 
the planet of 170°K. Since the scale height is to be used over the 
isothermal portion of the atmosphere, this author prefers to use the 
isothermal temperature (86 °K) . This drops the figure down consider- 
ably, and in fact makes i for model b = 26,000 feet. This compares 

with the value for Earth of 23,500 feet. 

The density profiles are plotted in Figure 2 based on the above 
assumptions. We should note that the density for model b is almost 
identical with Earth. Therefore, the dynamic pressures are equal, 
and ary winged vehicle which can fly above the Earth can also fly 
above Jupiter. Of course the pressures and lift required at a 
given altitude will be different. 

1 . 6 Summary 

Fairly conclusive evidence has been offered to support model b 
over mi>del a. Comparing the surface of Earth to the tops of the 

‘'•Sin 

clouds for Jupiter, one finds the density profiles nearly identical. 
Assuming an entry density at 160 kilometers for Earth, the same 
altitude would prevail for Jupiter. Even with model a, entry would 




(1-7) 



Here the scale height, ^ , is given by 




(1-8) 
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commence at about 25O kilometers. It is unlikely that any 
reasonable shape would penetrate below the clouds’ before the high 
Mach phase is over. Pressures increase very rapidly below the 
clouds, and the critical pressure for H2 is reached at about 
-120 km. From this point on we have solid to a pressure of 
200,000 atm,, and metalic beyond. It seems unlikely that 
Earth people would ever explore Jupiter to a significant 
depth below the clouds. But it is quite likely that they will 
explore the upper atmosphere for scientific reasons. Jupiter is 
the closest and easiest to reach of the remaining protoplanets. 
The Terrestrial planets all have secondary atmospheres, but an 
accurate determination of Jupiter’s atmosphere would give good 
evidence about the formation of the solar system. 

Such an entry could be made with a rather unsophisticated 
probe at an early date. The potential required would be very 
little more than the recently launched Venus probe. Higher g 
forces and heating would be encountered due only to the higher 
entry velocities involved. Such an entry will be designed in 
the subsequent chapters. 
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TABLE IV 

Summary of Atmospheric Characteristics 


Property 


Jupiter a 


Jupiter b 


Earth 


T 

o 


165°K 


168 °K 


288 °K 


P (slugs/ft. 3 ) 
* 0 


8.2 x 10’ 3 


-4 

9.02 x 10 


2.7 x 10" 3 


^(atm) 


24 


2 


1 


K (ft" 1 ) 


2.92 x 10"- 5 


3.84 x 10" 5 


4.26 x’lO" 5 


K (km" 1 ) 


.0958 


.1261 


>1397 

w. 


1/K (ft) 


34,300 

.#*' 


26,000 


23.500 


P(oK/km) 


2.6 


4.0 


6.5 


p. r (slugs /ft 3 ) 


1.915 X 10 " 3 


2.73 x 10" 4 


7.07 x 10" 4 


H tr (km) 


31.5 


21.0 


11.0 


H entry (kru) 


250 


165 


161 


gentry 


3 x 10 -12 


3 x 10" 12 


3 x 10” 12 


^iso 


86 °K 


86 °K 


217 °K 



I 
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CHAPTER 2 

GENERAL ENTRY DYNAMICS 



2.1 Introduction 

Since the purpose ‘‘of this thesis is to investigate design 
parameters for an unsophisticated probe to Jupiter for atmospheric 
determination, only trajectories which lend themselves to little 
or no guidance will be considered. This will limit us to a steep 
ballistic entry from deep space, a shallow entry from a disturbed 
reconnaissance orbit, or a decaying orbital entry. The use of a 
lifting vehicle will be prohibited by the delicate guidance 
required for orientation and attitude control. This first probe 
should be more like a meteor or bomb, and may broadcast only two 
things; stagnation point temperature and specific force. It may 
have to withstand specific forces up to 9000 Earth "g's", and will 
have to be dense enough to make a deep penetration before burning . 
It will use a thick ablative heat shield, and may rotate slowly to 
distribute heat evenly. It will contain no gyros or guidance 
system, and its transmitter will be potted solid state electronics 
and batteries. 

In all these one pass ballistic entries, the specific forces 
(deceleration) will be far greater than the gravitational forces 
on the vehicle, and therefore the entry can be considered a 
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straight line. A further assumption is made of planar motion. 

Since Jupiter rotates so fast, the orientation of this plane will 
be important. However, it should be relatively easy to make the 
entry near the equatorial plane, and in the direction of rotation. 

This will considerably reduce the entry velocity with respect to 
the atmosphere. 

The general three dimensional equations of motion for variable 
mass; oblate, rotating planet; and atmosphere with winds are 
spelled out in Duncan's thesis (ref 6). These equations cannot be 
solved in closed form, and are too complicated to be useful in 
showing variations of the parameters in question. For a probe as 
straight-forward as this, the two dimensional, constant mass, 
spherical rotating planet without winds will be entirely adequate. 

The vehicle is essentially a cannonball, and contains no guidance 
to utilize such nice refinements in motion. 

The ballistics used will be those of Chapman (ref 7) » and the 

modifications of Duncan (ref 6). Duncan's "conservation parameter" 

will not be necessary, since only the gas-dynamic phase is considered. 

The Keplerian and intermediate phases are of use only for guidance 

-4 -2 

in the range of 10 to 10 "g's", and do not affect vehicle design 

for the final phase. 

2.2 Types of Entry 

The direct entry will be made at some entry angle , , from deep 

space with initial velocity, v^ = v e sc* For ^ u Pi^ er this would be 
about 200,000 ft /sec. But if the entry were made in the direction 
of planetary rotation, this value would be reduced by 41,000 
cos 'y ft/sec. This type entry requires no rocket thrust, and no 
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guidance other than the initial aiming. 

The shallow entry would be accomplished by an impulsive 
perturbation from a stable reconnaissance orbit around Jupiter. 

The angle would be very shallow because large impulsive pertur- 
bations would require too much propellant. A reasonable value 
for ^ would be 5°» The entry velocity would be on the order of 
satellite velocity, and could be reduced by almost the entire 
rotational velocity of the planet. Therefore = 100,000 ft. /sec. 
This entiy would be gentler than the first type, but would require 
a higher ideal velocity capability for the launching vehicle, and 
much more accurate guidance. Extra fuel would be required to 
establish the orbit, and still more to launch the probe from the 
orbit . 

The orbital decay entry would be the gentlest possible with- 
out using a lifting vehicle, but it is also the most difficult 
to obtain. A low orbit within the sensible atmosphere would have 
to be attained, and would require the most careful guidance. The 
time required for this entry is longer than the others, as several 
passes around the planet would be made. The entry speed would also 
be about 100,000 ft/sec. and the maximum specific force would be 
40 Earth g's. This shows that a lifting vehicle of some sort will 
certainly be required if a manned entry is ever made. 

At the present time, the only entry available is the direct 
ballistic entry. But as we shall show, this would not be too 
hard to design with present materials, and would probably give the 
most useful information. 

Equations will be included for instrumentation of the probe 
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in both flight path axes and energy axes. Flight path axes can be 
used regardless of the shape or orientation of the vehicle. The 
assumption has been made that for a steep entry the specific force 
lies along the flight path (Drag). Therefore, pendulous mounted 
instruments will indicate flight path axes even in a spherical 
rotating vehicle. Energy axes (r, 0) can only be used, however, in 
conjunction with a vertical indicating system. The resolution of 
the motion into radial and tangential components requires the angle 
^ . It is expected that such a system will not be included due to 
its complexity and delicate nature. 

It should be emphasized that this instrumentation is to deter- 
mine deceleration forces for radio relay, and is not intended for 
guidance. In this matter, small angle approximation and neglect of 
gravitational terms is allowed. 

2.3 The Coordinate Frames 

Figure 3 illustrates the two simple two dimensional coordinate 
frames used in the analysis. The r, 0, { f' frame is a great circle, 
planet centered frame which rotates with the vehicle. It can be 
considered an air mass frame or energy frame, as potential energy 
and angular momentum can be written in terms of r and 0. The entry 
trajectory is assumed to lie entirely in the r,0 plane; and the angle 
of this plane with the equator, V' , is assumed to be small. 

The flight path coordinates (x,z) are in the plane of r,0 and 
are rotated through the angle, . 

2.4 The General Entry Equations 
Accelerations : 



A = R - 



r 



R 




( 2 - 1 ) 
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Components of Vj^ : 

Vj0 = V0 + Wjj R cos f = V cos ^ + Wjj R cos T (2-2) 
where v = vel. with respect to the air mass, 

and = inclination of the trajectory plane with respect to the 
equatorial plane. 

R = V sin (2-3) 

It is useful also to define coordinates along and normal to 
the velocity vector as an x,z frame. (See Fig 3«) The x„z frame 
is rotated from the r»jk frame by the angle ^ „ 

In x,z frame ? 

A x = V - W 2 -j-j R cos 2 T sin (2-4) 

A z = y ^ - V C ^ S — - W 2 ^ R cos 2 V cos - 2W JJ V cos V 



The forces to be considered are the gravitational forces and 
the drag force. 



G = - G 1 
sp r 

= - G (sin ^ 1 + cos 1 ) 

sp b x z 



(2-5) 



and drag ; 




Therefore the two dimensional equations of motion in the x,z plane 
are: 

V + (G sp - W 2 ^ R cos ^ ^ ^ sin + S = 0 force ) 

(2-7) 

vir- ^ c ° s b + (G - W 2 - T R cos 2 {^) cos - 2 VW TT cos 5/^=0 
b R sp IJ ' b IJ 

(z force) 
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Now, in dimensionless form, normalizing with respect to Rj, V 

I * si t • 

R j * ' 

tL+ - rA 2 cos 2 ^) sin + N D = 0 
r 

2 

v - 2L_ cos ^ + (—• - rA- 2 cos 2 5^) cos ^ 

In r, 0 coordinates: 
dv v 2 

-7-7= + N sin + (~ - r/L 2 cos 2 !^) - 2v<* ilcos Y - 0 

d T r D b J2. V 

r (2-9) 

+ + n d cos ^ + 2v r A cos ^ = 0 



( 2 - 8 ) 

- 2v/lcos Y - 0 



These equations can also 
angular momentum, rather than 
potential energy with respect 
the energy at infinity, total 

£l£eH = _ g r + 

M sp 2 



be written in terms of energy and 
coordinate forces. Defining zero 
to Jupiter's gravitational field as 
energy will be negative. 

( 2 - 10 ) 



In dimensionless form: 

v 2 

E = £(tot) = _I _ 1 

(tot) wy-2 2 “ r 

sat 

dE 

Equating to power : 



dE j 
dT" 



v. 



( 2 - 11 ) 



( 2 - 12 ) 



where f = the sum of all external forces in mean surface g's of 
Jupiter (dimensionless). This means of expressing motion is best 
for instrumentation, because accelerometers mounted within the body 
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can then measure the aerodynamic forces. 

HTT — 

“ = Vj (the output of the accelerometer in the 1 

direction) + Vj (the output of the accelerometer in the 
10 direction). (2-13) 

If the body has aerodynamic stability, more useful information 
would come from body fixed axes. 

£jS.= Vj (output of accelerometer parallel to V) 

X 

+ v^ (output of accelerometer normal to V) ( 2 - 14 ) 

The angular momentum of the vehicle about Jupiter provides the 
remaining energy expression needed. 

P = RV 

H 

non-dimensionalizing: p = 5—^ 

J sat 

= r (v^ + r/lcos V) (2-15) 

= F x f = r (accelerometer output in direction lx) (2-16) 
dT v 

If r can be measured by radar or integrated from other data, the 
equations of motion can be written in energy and angular momentum 
form; 

2 

dE _ d_ t v I 1 ) 
dT dT 2 ~ r 

= (Np)(v'+ rjicos ^ cos VO (energy) 

- 2 . = - rN^ cos ^ (angular momentum) ( 2 - 17 ) 

However, to utilize this system of instrumentation, the vertical 
will have to be known to provide the angle, <%£ ; and to establish 
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the r,0 coord, system model. We note that it need not be an accurate 
vertical as would be necessary for inertial navigation, but only an 
approximation to the angle, • This could be obtained from an 
"albedo tracker", or perhaps some other form of optical or magnetic 
device. 

Since refined measurements of forces are not needed, it is 
anticipated that the following data will be radioed back to Earth: 

1) The drag load factor, , which is the output of a 
pendulously mounted accelerometer within the vehicle. 

2) The surface or internal temperature, or if the vehicle 
does not rotate, the stagnation point temperature. 

2.5 General Solutions for Ballistic Entry 

These solutions are all in non-dimensional parameters as 
defined in Appendix I.. They are based on Duncan's approximations 
to Chapman's Z function analysis. Limits to the validity of each 
equation are stated with the equation. 

The basic assumption in the ballistic entry is that ^ is 
constant. That is, the., drag specific forces are much greater than 
the gravitational forces. This assumption is valid for velocities 
above 25,000 ft /sec in the vicinity of Jupiter, and for ^ such that: 

k^V >2.5 (2-18) 

o 

For Jupiter, this means an angle- greater than approximately 2.3°. 



dh 

dT 



- v sin 



(2-19) 



d<r 

dT 



- tocsin r, 

dh dT o 



( 2 - 20 ) 



The non-dimensional drag load factor, N^, is defined as: 
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c s 



■» -- fr- I /*• R , -5- * 



J M 






( 2 - 21 ) 



Therefore, 

dv = dv dT _ 1 /°o R 

do- ~ d T dcr '2k sin ^ M 



V 



( 2 - 22 ) 



Separating variables and integrating, we have velocity as a 
function of density ratio, and atmospheric and vehicle parameters; 



v 1 /^oR C D S 
in ” — 



2 k sin^ M (O' - cn ) 



(2-23) 



Range angle in the ballistic entry is found by the approximations 



\ - \ = cot $ ln t . 4 oot 3s (h - h i ) 



(2-24) 



It has been assumed in the first chapter that: 
-kh 

<T = e 



(2-25) 



All these solutions for altitude, density, and velocity can be 
written in terms of each other, or in terms of dimensionless time, 
To 

Density ratio as the independent variables 



1) h(<r) = - Or) lncr 



2) v(tf-) = v i e 



2k sin^M 



Altitude as the independent variable 
1) <T(h) = e- 1 * 

°o R "D’ 



C n S 

D )(e 



2) v(h) = v,e 2k sin^TM 

p 



kh = e -**i) 



Velocity as the independent variables 



( 2 - 26 ) 

(2-27) 

(2-25) 

(2-28) 



(2-29) 
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2) h(v) 



1 , f2k sin^M „ / v \ -kh^ 

* - k ln |/>o R C D S 111 + e 1 



(2-30) 



Dimensionless time as the independent variable: 

1) <T(T): 



-kv. sin 
l c> 



.AsVa 

“'2K sin tf.vO 

p 



(T- rj 



= ln 



_.i co o R C_S n 

—I + e (_ L — _2_) n (£: — i.) 



<3T| n _^ ' 2k sin^M 



n • n i 



2) h (T): 



Po R CjjS 

-Vo 1 ,; e i ) 



(2-31) 



- tr-Ti)^ sin 2k sin ^ M 

= k (h, - h) + I (- aj^) (e-^ - e^i) 
n=l £ 



n • n 1 



(2-32) 



3) v (T): 



^oRC d S 

- (T- TJ kv ± sing'e 2k sin | M 



= ln 



2k sin$M , /V \ , 

7°ic^Kp 1 ln ( ~ > + 1 



V. 

1 



7 <- i" Irf 



n=l 



n 



n 



(2-33) 



The summation terms in the functions are the remaining terms from 
the series approximation to the exponentials , and may be dropped at 
the user’s discretion. 

We note that the vehicle parameter, M/C^S , turns up in every 
equation, as well as the atmospheric parameter , /^oRj /2k, and the 
entry angle . Specific solutions of the equations of motion for 
an entry into Jupiter are presented in the next chapter. These 
solutions will show the variation of the design parameters for 
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specific atmospheric parameters and entry angles, 

2,6 Solutions for Specific Force and Stagnation Temperature 
The two greatest problems in planetary entry are the high 
specific forces and heating encountered. 

The specific force for a lifting vehicle at high Mach is 
expressed in non-dimensional form ass 



. G J 2 A,R C D S 
E ~ G e V 2M 



O' 



, /Lx 2 
1 + (n) 



j 



1 1/2 



( 2 - 3*0 



This, in turn, can be expressed as a function of ary other variable 
by using the equations in section 2,5« For example, as a function 
of velocity for a non-lifting bodys 



„ , x g j 2 Po RC n s 

M * 5- v — ar 



zm j. 2k sin^M — 
i /O r r . _.<? n v. 



/’o RC D S 



(2-35) 



This expression is in Earth ”g 5 s", and is valid when << <J~» 
This function will have a maximum whens 



h (vV) ■ 0 



thus 

^.= > 2 ! sir. r, 

dT 2 0 

and 

v (f max,) = ,607 v„ (2-36) 

E X 

This result is independent of vehicle parameters. The 
expression for maximum specific force depends on atmospheric 
and entry parameters only. 

f^ (max) = (max) = — i 1— sin ^ 



(2-37) 
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The height at which this maximum occurs, however, is a function 
of the vehicle design parameter, M/C^S . 

There are several heating factors of interest, but the stag- 
nation point temperature will be the hottest to contend with. 

Also, the velocity at which the maximum stagnation point temperature 
is reached is independent of atmosphere and vehicle parameters. 



and "heating function" is an atmospheric constant relating Jupiter 
to the Earth as follows: 

ratios of: 



The value of this heating function for Jupiter is about 64.2. This 
will be discussed more fully in the next chapter. 

The stagnation temperature as a function of velocity is: 



T 



s 



= 1.392 x 10 4 (HF)£/ 4 (VF) 1 / 4 (vV) 1 / 8 



( 2 - 38 ) 



where "vehicle function" (VF) = p=- 

^rad v 



(2-39) 





(2-41) 



This is a maximum when <T~) = 0 

dv kv^ sin 

d'f Z 



and v(T max.) = .846 v 
s 



i 



(2-42) 
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Use will be made of these factors, to show where maximum 
specific force and tempferature are expected, in the velocity vs. 
altitude plots in the next chapter. 

Little has been said of the decaying orbital trajectory. 

Such an entry would require an extreme accuracy in guidance, 
beyond the scope of this simple probe. The escape from Earth 
and establishment of a reasonable transfer ellipse (2 yr 9 mo.) 
requires an ideal velocity potential of about 14.3 km/sec (ref 5)» 
To establish a stable orbit about Jupiter at all will require an 
additional ideal velocity of 5»7 km/sec. This would be taxing 
the present level of rocketry to its limit. A totally rocket 
braked soft landing would require 66 km/sec. more, and is out 
of the question. 

The equations for establishing a decaying orbit are out- 
lined in Duncan (ref 6). There is first the circularization 
phase during partial passes through the atmosphere; then the 
decay phase. These equations are non-linear, and have not been 
solved in closed form. However, machine solutions give a maximum 
Np for Jupiter of 40 Earth g”s, and a stagnation point temperature 
of 2.7 times that for entry into Earth (ref 7) • It appears that 
the unmanned entry is possible, but the establishment of the 
necessary low orbit is beyond the state of the art in guidance. 

We will therefore concentrate on the unguided ballistic entry. 
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CHAPTER 3 

BALLISTIC ENTRY INTO JUPITER 



3.1 Introduction 

This chapter deals with the specific entry profiles into 
Jupiter. Solutions to the equations presented in the second chapter 
are calculated for entry into the two model atmospheres. For each 
atmosphere, and for several entry angles and velocities, velocity 
profiles are calculated. These show a variation with the vehicle 
design parameter, M/C^S. 

Also developed here are specific force and heating profiles. 

The heating problem will be considered in two phases : 

1) Total heat absorption, Q^.. 

O 

2) Maximum heating rate, • 

Preliminary calculations have shown a range in the vehicle 
design parameter from 4 to 100 slugs/ft. , and these values are 
included in the figures. Finally, with an eye toward instrumentation, 
specific force and stagnation heating rate plots against real time 
are generated. These vary with design parameter, atmosphere model, 
and entry conditions, and are the predictions upon which experi- 
mental atmospheric determination is suggested. 

3.2 Velocity vs. Altitude 



The expression as developed in chapter two is calculated here: 
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.1* 

v (h) = 2k sin <£ M * ' 



-kh. \ 
e l) 



(3-D 



The plots are made as a ratio of v/v^ against altitude, H. There- 
fore, an entry from any initial velocity can be considered from 
these graphs. The initial conditions vary with the model atmos- 
phere as follows ; 

1) Model a. Figures 5 a * c 

a) = ,0082 slugs/ft. 3 

b) k = KR 

J * 

= 7615 

c) M/CpS = 4 to 100 slugs /ft, ^ 

d) (^ = -5°, -30°, -60° 

e) h. = H./R, 

i l J 

= ,003 (H t = 500,000 ft.) 

2) Model b. Figures 4a, b, c 

a) /° = .0009 slugs /ft. 3 

b) k = 8840 

c) M/C d S = 4 to 100 slugs /ft. 2 

d) ^ = -5°, -30°, -60° 

e) h ± = .002 (H = 333,000 ft.) 

The values of ^ were chosen as representative of the direct 
entry, and also the disturbed orbital entry. As previously 
stated, an economical use of mass ratio will only allow the 

= -5° case for disturbed orbital entry. These solutions are 
valid to about .25 v^. Beyond this point, gravitational forces 
are significant, and changes in Cp and other parameters are 
apparent as Mach number falls below about ten. However, the 
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profiles cover the range in which heating and deceleration effects 
are greatest, and will therefore be useful. These solutions were 
compared to sample calculations by machine from Chapman (ref 7)* 
and agree closely with the more accurate solutions obtained there. 

The most significant finding here is that an extremely dense 
vehicle will be required to reach the cloud layer, if model "a" is 
accepted. Even at <3^ = -60°, M/C^S must be 100 to arrive at the 
cloud layer while still in the entry phase. If we wish to ascertain 
the properties of this cloud deck, we must arrive there while there 
is still a clean cut response from either the heating rate or 
specific force measurements. This is easily accomplished in the 
model "b" atmosphere. 

3.3 Specific Force 

The specific force in earth g's has been given in chapter two as 




(<T + 2k sin <^M 
1 / ° 0 ] * C D S 




(2-35) 



This force is assumed to be all along the flight path over the 
portion of the trajectory where G is not considered. In addition, 
if we consider v i as entry velocity, and <71 as negligible, the 
maximum specific force is independent of vehicle design . Thus the 
non-dimensional drag load factor is given by: 

Q 

N d (v) = k sin <4^ v 2 In 2- (3-2) 

E i 



This load factor is a function of entry angle, entry velocity, 
and present velocity. Present velocity, as a function of height, is*- 
also a function of vehicle design. The maximum load factor, however. 



is a function only of the atmosphere, entry angle, and entry velocity 
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The entry velocity relative to the air mass is in turn a 
function of entry angle and planetary rotation. In the case of 
the Terrestrial planets, the rotation is negligible. But Jupiter's 
large radius and fast rotation period (9 hours, 5 0 minutes), make 
inclusion of rotation velocity necessary. Figure 6 shows a- plot 
for both atmospheres of vs entry angle considering rotation for 
the direct entry. Jupiter's equator lies nearly in the ecliptic 
plane, and we assume a probe will remain in the ecliptic plane for 

minimum energy requirements. Therefore, a range of 180° in entry 

•* 

angle along the equator is covered, considering the preferred 
entry with the rotation, and the disastrous one against the rota- 
tion. The graph clearly shows the effect of the rotational 
velocity of 41,000 ft/sec. 

We note with some concern that for the usual entry, atmosphere 
b with the rotation, we can expect 100 earth g's per degree of 
entry angle'. Specially constructed vacuum tubes can take about 
400 g's. Therefore we must limit the probe to 4° entry angle if 
we use a vacuum tube transmitter. This limit is too restrictive 
on the guidance system, and we will have to use solid state 
electronics capable of about 3000 g's for a 30-degree entry. 

Accelerometers are available to measure up to 10,000 g's 
(piezo-electric accelerometers). The accelerometer system and its 
two gimbal mount with four roller bearings will be designed for 
10,000 g's. We note for interest that for a .1 pound accelerometer, 
the outer gimbal bearings will be freight car type roller bearings 
with 50 , 000 pound capacity each. 

For a disturbed orbital entry, the maximum expected is 
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shown in Figure ?. Here the entry angle is more easily controlled, 
and the entry velocity is lower. The expected is 40 g's per 
degree entry angle. Also, because of fuel limitations, a maximum 
entry angle of five to seven degrees is expected. The vehicle could 
then use vacuum tube electronics, and be designed for 250 g's. But 
all this relief comes at great expense when it is realized that the 
vehicle must first be put into a stable orbit about Jupiter. Also, 
we see from recent experience with the Discoverer satellite series, 

i 

that firing a retro-rocket by remote control is largely unsuccessful. 

The variation in specific force with altitude was calculated, 
and plotted in Figure 8. This graph shows how the vehicle design 
parameter controls the height for maximum N^. There are plots for 
the several entry angles used as examples for parabolic and orbital 
decay entries. This information can be used to determine the 
vs time profile for the selected design. It represents dv/dt vs H, 
and therefore H vs time can be derived. It is rather too much work 
to derive this for every possible case, but it will be done for the 
design chosen in the next chapter. The specific force variation 
with time will be one of the transmitted outputs. 

Figure 9 is an illustration of the time profiles to be expected. 

All quantities are specified, and the results in N^, H, and V are 
plotted. We see here that the entry phase is veiy short compared to 
an entry into Earth. The deceleration phase takes approximately 
three minutes for Earth, while it only lasts 15 seconds for Jupiter. 

This short time will aid us in the structural considerations for design. 
3.4 Heating 

The expected heat generation profiles will be calculated here. 
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The problem of handling the heat will be considered with the 
vehicle design. Most sources agree that for a ballistic entry, 
re-radiation of heat is insignificant. These people are thinking 
in terms of an Earth entry of about three minutes duration. For 
the Jupiter entry of 15 seconds duration, there will be even less 
chance of a radiative equilibrium. We will, therefore, only 
consider the convective terms in the heat flow equations as 
presented by Chapman and Duncan. 

This convective heat will occur almost entirely in the 
continuum flow regime. According to the Blasius boundary layer 
equation, the heat transfer through a laminar boundary layer is 
inversely proportional to the square root of the local Reynolds 
number based on the radius of curvature of the nose. Thus, for 
bodies in the same environment , the heat input rate per unit area 
and time will be inversely proportional to the square root of the 
radius curvature of the missile nose. 

In the ablative heat shield case, a sharp nosed body would 
soon be reduced to a blunt nosed body due to the small amount of 
ablative material in the pointed portion of the missile. The 
primary reason for using a bLunt body is to secure the high Cp 
necessary for deceleration at high altitude. 

The stagnation point convective heat rate is given by the 
expression: 



O' atm is a function of the atmosphere, and is given as 18,000 for 
Earth, and works out to 1,590,000 for model b and 4,750 ,000 for 



(Q ) 
v c s 




C atm ,/<>s n , V ^ BTU 



( 3 - 3 ) 
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model a. Compared to Earth, the heating rate will be 87 and 264 
times worse, respectively. 

In terms of our non-dimensional parameters, the heating rate: 

Q = «T} n (V) m HS (3-4) 

s VrcT ft 2 -sec 

a 

The exponents on density ratio and velocity ratio depend on 
the flow conditions in the boundary layer. This in turn depends on 
the Reynolds number. 

For laminar flow, where p, is proportional to V T, then n = l/2 
and m = 3* 

For turbulent flow, n = .8 and m = 2.2 (from Chapman). Figure 
10 is a plot of Reynolds number at maximum stagnation temperature 
velocity (.846 v ) for various M/C^S ratios. Chapman states that 
we can have high hopes' of maintaining laminar flow in the boundary 
layer. But this plot shows that laminar flow can hardly be expected 
over much more than one foot of the Jupiter probe vehicle. If we 
assume the boundary layer will become turbulent at 6 or 7 million, 
a body of M/CpS = 10 at an entry angle of 30 degrees will have 
laminar flow for about one foot from its leading edge. However, 
if the body is blunt enough (e.g., like a Mercury Capsule) most of 
the heating will occur over the laminar flow regime. 

These extremely blunt bodies will suffer more from radiation 
heating from the hot gas to the vehicle, but the time of flight is 
so short we believe radiation will have little effect on the total 
heating problem. 

The maximum laminar heating rate is given by the expression: 

<Vs <■*“> <- 605 > (cr)1/2 ‘V 3 



(3-5) 
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This is a function of M/C^S and entry conditions. The density 
ratio used is that for which v = .846 v^, and is calculated from 
the graphs in Figures 4 and 5° Figure U is a plot of the maximum 
laminar stagnation point heating rate expected vs. the design 
parameter M/C^S. Also shown on this graph for comparison are the 
maximum allowable heating rates for Teflon and graphite ablative 
heat shields on a spherical vehicle (ref 8,9). This graph 
indicates at once that we will not be abLe to tolerate steep 
ballistic entries with our present ablative protection. This is 
one of the important constraints to be used in the next chapter. 

Figure 12 is an example of a heat rate vs. time plot for a 
particular entry configuration. The integral of plots such as 
these will give the total heat absorbed, and will determine the 
weight of ablative material to be used in the vehicle design. 

3.5 Summary 

We have now developed particular constraints for the Jupiter 
probe. The structural limitations under high specific force, and 
the heat dissipation problem will determine the upper limit of 
severity of the entry. Fortunately, human factors are not involved. 
The increase in guidance requirements, and mass ratio of the 
booster will determine the practical floor on the gentleness of 
the entry. In the next chapter, we shall see whether these 
constraints leave a channel of possibility for the Jupiter probe 
using present day hardware. 

The heating rates developed here are for a laminar boundary 
layer. The heating rates in a turbulent boundary layer, due to 
better convection, would be prohibitive. But we must not fool 
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ourselves into thinking the boundary layer will remain laminar very 
far aft of the nose. This requirement on boundary layer will present 
a further constraint. 

Stagnation point temperature, per se, has little relation to 
the entry profile when an ablative coating is used. The shield 
quickly heats to the vaporization temperature, and remains constant 
while the heat is dissipated. This boiling off of the gaseous 
material further shields the vehicle itself from convective heat 
input. The vehicle skin will remain below the critical 2000°F 
temperature (for steel) as long as the coating lasts. This is due 
primarily to the poor heat conduction of the ablative coating. 

These considerations will influence the instrumentation of the probe. 

The problem of communication of the data back to Earth is not 
considered in this thesis. However, it should be recognized that an 
ionized shock layer may interfere with electromagnetic transmission 
at certain frequencies. The detail design of the antenna, and the 
selection of transmitting frequency, will have to include considera- 
tion of the degree of ionization in the atmosphere around the vehicle. 
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CHAPTER 4 

CONSTRAINTS ON VEHICLE DESIGN AND FEASIBILITY 
OF ATMOSPHERIC DETERMINATION 



4 . 1 Introduc ti on 

At this point, we must examine the constraints imposed on the 
vehicle design by the forces and heating rates calculated in the 
previous chapters. We will then estimate the constraints due to 
the present limitations in guidance. The vehicle configuration 
will be specified. Within the area of possibility for a success- 
ful entry, we will discuss the feasibility of atmospheric model 
determination. The conclusions and recommendations will follow. 

4.2 Heating Rate Constraints 

We have discovered that the only feasible means of heat 
dissipation will be an ablative shield. In view of this, and the 
short time duration of the heat input (2 to 10 seconds), the heating 
of the vehicle core should be acceptable. We can carry enough 
heat shield without a great payload penalty to handle even this 
large amount of heat. The vaporization of the material will not 
only carry away heat to the surrounding air, but will provide an 
insulating region of vaporized material to further shield the 
vehicle from heat convection. Despite the large heat inputs, the 
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materials used are poor conductors, and the interior should remain 
below its critical strength temperature (2000°F for steel) . 

The material most generally used for this type shield is Teflon 
(a Dupont trademark for polyethoflourethylene) . It is a plastic 
with strength at high temperatures, and is easily molded. Its heat 
capacity to vaporize is about 10,000 BTU/lb., and the maximum heat 

r\ 

rate input is about 50*000 BTU /ft. sec. (ref. 8). This limit keeps 
the material from boiling off so fast that it explodes. 

The best material for an ablative heat shield from the point of 
view of heat capacity is graphite. Graphite can dissipate about 
30,000 BTU/lbs., and at a rate of 150,000 BTU/ft.^sec. But there 
are many structural problems associated with the use of graphite 
that have not yet been solved. It is brittle, and hard to cast 
unless it is mixed with other more plastic compounds. But such a 
mixture lowers the heat capacity. Perhaps a mixture of Teflon and 
graphite may be successfully used. 

The heat rate input graph (Figure 11) shows the limits imposed 
by these materials on the type of possible entry. The Teflon 
coated sphere of reasonable density (from 1 to 4 with respect to 
water) can only be used for the deflected orbital entry. If we 
allow a solution to the graphite structural problem, we can still 
only make a '20 degree direct entry. These are the rather stringent 
constraints imposed by the maximum allowable heat rates. 

4.3 Configuration 

Another glance at Figure 11 shows us that the lighter vehicles 
will allow a steeper trajectory than the heavy ones. Most of the 
missiles used today have a payload specific weight very close to 
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one (with respect to water) . Figure 13 shows the variation in 
M/C d S with radius of curvature for a sphere. Figure 14 shows this 
information for a typical final booster stage, and Figure 15, for 
a project Mercury shaped capsule. A of 1.0 is used for these 
shapes at hypersonic velocities. The Pioneer V probe has shown a 
deep space launching capability of about 100 pounds payload. The 
Saturn booster is reported to have a 15,000 pound payload. These 
lines are also indicated on the configuration graphs. 

It is apparent that if we limit ourselves to 100 pounds , we 

r\ 

can only attain M/C^S of about 4 slugs /ft. . This would be to our 
advantage for heat protection, but would not allow penetration to 
the vicinity of the cloud level during the entry phase. Using the 
Saturn booster, we can attain M/C^S on the order of 20 slugs /ft. ^ 
This would give us the desired profile in atmosphere model b. We 
would have to go to still greater weights to reach the clouds for 
atmosphere model a. 

We note that the project Mercury qapsule has the advantage of 
a larger radius of curvature for its diameter than the other 
typical shapes considered. This shape has been demonstrated to 
have aerodynamic stability at high Mach number, but a still flatter 
nose would probably lead to instability. In Figure 15 we have run 
the calculations to a density of 7» This would be a solid steel 
vehicle with a graphite heat shield, and would approximate the 
limit of density in vehicle design. 

In order to show variation of design to best advantage in 
terms of mission, entry angle, and M/C^S we must assume some 
basic limits. We assume that the Saturn booster is available. 
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and that a graphite heat shield can be used. We therefore let the 
weight equal 15 , 000 pounds. We allow average vehicle density to run 
from 7 to 0.3 with respect to water. We limit the vehicle size to 
a 10-foot diameter. 

Using these criteria, and satisfying the maximum heat rate 
constraint, we develop the plot in Figure 16. The variation of M/C^S 
with entry angle as a function of vehicle average density for the 
three shapes chosen is shown. The mission requirement is arbitrarily 
set as reaching the clouds with v = .25 v^. From Figure 1 6, several 
facts are evident. 

1) With the very low M/C^S compatible with the present 100 
pound launch capability, the vehicle cannot hope to reach low 
altitudes during the entry phase. We must therefore wait for the 
Saturn booster. 

2) The cortifcal shaped vehicle has an advantage due to its 
larger radius of curvature for a given vehicle diameter. This 
vehicle can come fairly close to meeting the mission requirement. 

3) It would tax the imagination to conceive of a vehicle 
large enough to penetrate atmosphere model a to the cloud level. 

r\ 

M/CpS would have to be 100 slugs/ft. , and with a solid steel vehicle, 
the diameter would be 20 feet, and the weight would be 127,000 pounds. 
In addition, the ablative material would have to absorb 450,000 
BTU/ft.^ sec. 

4) There will be a lower limit on the vehicle average density 
dictated by the specific force requirements. A detailed structural 
analysis is beyond the scope of this thesis. The specific force 
goes up practically linearly with entry angle if we enter in the 
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direction of planetary rotation, and figures to about 100 Earth 
g's per degree. We suggest that a reasonable structural limit for 
a design load factor of 3000 would produce an average vehicle 
density of two. 

5) It is evident from the figure that high M/C^S alone does 
not guarantee a deep penetration. The heating rate constraint 
limits such dense vehicles to a low entry angle for parabolic 
entry velocity. 

6) The deflected orbital entry is virtually free from heating 
constraints. Any of these vehicles with an M/C^S of 100' can enter 
at more than 30° without overheating. But the deflected orbital 
entry is not that easily accomplished. We are limited by fuel 
carried in the vehicle to a deflection of 5 or 10 degrees. From 
Figure 16, we must have a high M/C^S. At 5° deflection we must 
have 100 slugs/ft. Only the rocket shaped vehicle can reach 

that high an M/C^S without exceeding the 15.000 pound weight limit. 

The configuration chosen as most likely to succeed is the 

Mercury capsule shape with an average density of two. The entry 

angle should be near 30° • The M/C^S is 11, the nose radius and 

vehicle diameter is 8.5 feet, and the weight is 15,000 pounds. 

A graphite heat shield is used, and the maximum heating rate w ill 

be 150,000 BTU/ft.^sec. The significant heating period is 2.4 

seconds, and by graphical integration the total heat absorbed is 
2 

180,000 BTU/ft. of heat shield (Figure 12). This specifies the 

4 

2 

shield weight as 6 pounds /ft. (Ref. 9), and makes total ablative 
material weight 1360 pounds, or about 9$ of the vehicle weight. 
Other materials or cooling systems would require much greater 
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weight, and this 9$ seems to be a practical low for the Jupiter entry 
vehicle. The interior will stay relatively cool as long as there 
is some of the ablative shield remaining. If we make the shield 
lighter, the structure would melt or lose strength as soon as the 
ablative material evaporates. The majority of the deceleration 
loads occur after the heating phase, and we must not allow the vehicle 
to become weakened early in the trajectory. 

The heating protection described here is for the model b 
atmosphere. From Figure 11 it is evident that only a 5° entry 
could possibly be tolerated in atmosphere a. Such an entry would not 
meet mission requirements, or guidance constraints. It is only 
practical, then, to » design the vehicle for model b, and allow its 
failure in model a to be an indication of the atmosphere. 

4.4 Guidance Constraints 

It is difficult to say just how good our present guidance 
systems are. The Pioneer V probe to the orbit of Venus was judged 
to have performed perfectly. Yet it is short of its goal by several 
thousand miles due to a booster burning time uncertainty of less than 
one second; It seems clear, then, for a precision injection at any 
controlled angle into Jupiter, we must have mid-course and terminal 
guidance. Such a system, furthermore, must operate automatically 
without signals from Earth. It must take its measurements in the 
vicinity of Jupiter from the fields of Jupiter (optical, magnetic, 
or electric). It would probably carry a model inertial reference 
frame (stable platform with gyros), which would be corrected for 
drift periodically by star trackers during its 2 year 9 month journey. 

This system does not exist today, as far as the unclassified 






literature is concerned. If it did, it would have been used in the 
Pioneer V probe to Venus. Such developments may reasonably be 
expected in the next ten years, however. The consensus at this 
Institute is that with our present launch guidance we would be 
lucky to hit Jupiter at all. The Keplerian trajectory and ballistic 
entry make no provision for corrections to launch errors, and 
certainly we could not expect a precise injection angle, much less 
a stable orbit about Jupiter. 

It might be interesting to consider the statistical chances 
of attaining a 30 ° or less injection in the direction of planetary 
rotation, assuming we can attain a parabolic contact with the 
planet. This would be an academic study at best, for we must not 
expect the United States to expend this awful amount of money 
without a good probability for success. It seems we are con- 
strained to wait for a system including terminal guidance. 

4.5 Feasibility of Atmospheric Determination 

For this discussion, we shall assume we can make the 30° 
entry with the vehicle selected under section 4.3. What informa- 
tions hould we instrument and radio back to Earth? 

An examination of Figure 6 indicates that there is considerable 
difference in the maximum specific force expected in the two 
atmospheres. Model a will present an (max) 2ty % lower than 
model b. This difference is quite sensitive to entry angle, how- 
ever, and an error of 10° in entry angle will also equal 24$ of 
Np. We must therefore include instrumentation to measure the 
entry angle. Even if this angle cannot be controlled, the model 
atmosphere can be ascertained if the angle and its degree of 



certainty are known. The maximum specific force is a definite 
quantity not dependent on the vehicle configuration, and therefore 
not subject to uncertainties such as C^, or effective drag area 
(which may be considerably larger than the cross sectional area of 
a vehicle in an ionized gas). This force is a function solely 
determined by the entry parameters (which we have assumed are 
successfully controlled), and the atmospheric parameter, k, to be 
determined . 

The skin temperature or its derivative appear to hold atmospheric 
information for us also. Figure 10 shows that the Reynolds numbers 
for the two models are nearly identical. However, the heating 
profiles are quite different. The higher cloud level pressure and 
density of model a account for the higher heating rates (Figure 11). 
Unfortunately, a temperature would be difficult to instrument in a 
vehicle incorporating an ablative heat shield. The heating is wholly 
transient, and in the short effective heating period, no equili- 
brium temperatures are reached. Further, the temperature or heating 
rate are quite sensitive to uncertainties such as C^, S, Prandtl 
number and viscosity. To use heating as a determinant, these 
instrumentation difficulties and uncertainties would have to be 
overcome . 

A very useful determinant would be altitude, and its measure- 
ment by radar would also provide the entry angle needed for 
determination. The height at which maximum occurs varies 
considerably with the atmospheric parameter (Figure 8) . The clouds 
of Jupiter are expected to be ammonia cirrus (solid crystals), and 
would give a good radar echo to the proper frequency emission. The 
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radar could be made to track minimum range, and the dish angle 
would then be the flight path angle. 

There would be difficulties with the radar, of course. The 
tracking mechanism would have to be perfectly balanced and very 
rugged to endure 3000 g's. The electronics would have to be all 
solid state. There may be interference, as with the radio trans- 
mission, from the ionized shock layer. But these problems are not 
insurmountable, and should require no more development time than 
the guidance system required. 

Even without the radar, specific force measurements would be 
indicative of the atmosphere to the degree of uncertainty of the 
entry angle. 

4.6 Conclusions and Recommendations 

The major problem appears to be control of the entry angle. 
Atmospheric determination can be made if the entry is successful, 
and if the entry angle is known within about five degrees. To 
make a successful entry, the entry angle must be controlled to 
30° or less. 

The mission is stymied at present by the lack of a mid- 
course and terminal guidance system that would provide this 
accuracy in entry angle. Such a system is probably under develop- 
ment at this time. It need not be stressed for the high load 
factors encountered during the entry phase. It would be sufficient 
to control the entry angle prior to the onset of the specific forces. 

If this controlled entry angle were reliable, the only 
instrumentation needed in the entry vehicle would be a pendulously 
mounted piezo-electric accelerometer, A potted electronics 
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transmitter could be built to broadcast the load factor thus 
measured. If the controlled entry were not reliable to within 
five degrees, a device would be needed in the probe to measure the 
actual angle. If radar is developed for this measurement, it can 
also supply the vehicle height above the clouds. This height can 
be used as a cross check in atmospheric determination. 

The vehicle itself will have to be large. In the Mercury 
capsule configuration ' the vehicle could reach nearly to the cloud 
level at a weight of 15,000 pounds. This will require the Saturn 
booster for launch and a graphite heat shield — both under develop- 
ment at this time. The present deep space launch capability of 100 
pounds will not be adequate for the Jupiter probe. A Teflon heat 
shield does not have enough heat rate capacity for the direct entry, 
but would do for the deflected orbital entry. 

The decaying orbital and deflected orbital entries are still 
farther into the future with respect to guidance. 

The primary purpose of the mission is to determine the molecular 
weight of Jupiter's atmosphere. It has been shown that this can 
be done with a large, but simple, entry vehicle. But the entry 
vehicle must be delivered to the entry point by a much more sophis- 
ticated navigational vehicle. This delivery is beyond the present 
state of the art of guidance. 

We suggest for further study the following points : 

1) The detailed structural design of the vehicle to withstand 
the high load factors encountered; 

2) The degree of ionization in the shock layer, and its effect 
on electro-magnetic transmission; 
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3) The chemical effect of the heat and the reducing atmos- 
phere on the ablating material; 

4) The probability of achieving a 30° or less entry angle 
with present day launch guidance only. 

We believe the guidance, booster, and graphite heat shield 
problems are currently under consideration. 
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A. Subscripts 

r. </>,r 

i 

I 

J 

E 

b 

s 

B. Symbols 
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^atm 

C D 

°P 

°V 

D 

s 

E 

F 

f 



G 

I 

H 



APPENDIX I 

NOTATION AND SYMBOLS 



Components along the ly , , 1 ^ axes. 

Initial value of the quantity. 

Measured with respect to inertial coordinates. 
The quantity for Jupiter. 

The quantity for Earth. 

Refers to ballistic entry. 

Refers to stagnation area. 



Acceleration vector of entry vehicle, ft/sec^. 

Dimensional constant for the planetary atmosphere. 

_ . (BTU ft" 3 / 2 sec" 1 ) 

Drag coefficient . 

Heat capacity at constant pressure. 

Heat capacity at constant volume. 

Drag force (lb.). 

Energy (ft-lb). 

J? 

Dimensionless energy = — “ — 

MV 2 

sat 

Specific force (ft/sec. ) 

Dimensionless specific force in surface g's of 



- F 

the planet. f = -g- 

2 

Gravitational field intensity (ft/sec ) 

2 

Planet's gravity field (ft /sec ) 

Altitude of vehicle above planetary reference level (ft). 
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h 

K 



K rad 

k 



M 



M, 






P 

P 



RC 

R 

g 

R 

T 



Re 

S 



T 

t 

V 



H 

Dimensionless altitude = 

K 

Exponential decay parameter of planetary atmosphere 
(ft -1 ). 

Vehicle surface radiation emissivity. 



Dimensionless decay parameter of planetary 

atmosphere = KR 

Mass of vehicle (slugs) 

Mass of planet (Jupiter, in this case) 
Dimensionless drag load factor = 

Angular momentum per unit mass of the vehicle 

p 

Dimensionless angular momentum = — rrjx 

R*'* O' ' 



Prandtl number 

Total convective heat absorbed per 'unit area. 




Radius of curvature of vehicle nose (ft). 

Universal gas constant = 8 . 31 x 10? ergs/°K 

Radius vector from planet center to vehicle 

Dimensionless radius vector = = — 

R j 

Reynolds number length 

M- 

Reference area of entry vehicle used in drag 
computations (ft ). 

T emper a tur e ( °K ) . 

Real time (sec). 

Velocity vector of a vehicle with respect to 
coordinates rotating with a planet. 
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Dimensionless velocity = ^ 



V 

sat 



V 



sat Circular orbital velocity at planet surface = -yGR 
X Distance flown measured at surface of planet 
Z Chapman's transformation variable = v cos £b s ^ n ^b 

k 



C . Greek Symbols 



y Ratio of specific heats Cp/C r 



y b 

r 

/’o 

<r 



Ballistic entry angle 

Dry adiabatic lapse rate (°K/km) 

Coefficient of viscosity of the planet's atmosphere 
( slugs /ft-sec) 

Free stream atmospheric density (slugs /ft'*) 
Reference level atmospheric density 
Density ratio /° I/O 



S( ) Summation of quantities in parenthesis 



n=l 

T 

0 

r 

n 



Dimensionless time = 



sat 

R 



Angle measured in plane of trajectory from a 
reference point in the direction of motion 
Angle of inclination of the trajectory plane with 
the planet's equatorial plane. 

Dimensionless angular velocity of the planet about 
its polar axis 



n - *J¥ 



for Jupiter 
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Xi initial range angle 



ir 




FIGURE 3 ENTRY COORDINATE SYSTEMS 
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